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Abstract 



We suggest that photons with energies between 5 and 10 MeV, generated by 
the {u, u'p'j) and (z/, v'n^) reactions on ^^O, constitute a signal which allows 
a unique identification of supernova and Vr neutrinos in water Cerenkov 
detectors. We calculate the yield of such 7 events and estimate that a few 
hundred of them would be detected in Superkamiokande for a supernova at 
10 kpc distance. 
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Neutrinos play a decisive role in various stages of supernova evolution In particular, 
the gravitational binding energy of the nascent neutron star is released by neutrino-pair pro- 
duction [0] . It is the neutrinos generated during this cooling and deleptonization phase of the 
hot remnant core which will be mainly observed in earthbound detectors. Although pairs of 
all three flavors are generated with equal luminosity due to their smaller opacities t'^ and 
Ut- neutrinos and their antiparticles decouple at smaller radii, and thus higher temperatures 
in the core, than Ue and Ue neutrinos. As the neutrinos decouple in neutron-rich matter, 
which is less transparent for than for z/g, it is expected on general grounds that the neu- 
trino spectra after decoupling obey the temperature hierarchy [^], T^^ > T^^ > T^^, where 
stands for z/^, 1/^- and their antiparticles, which are assumed to have identical spectra. The 
neutrino spectra can be approximately described by Fermi-Dirac (FD) distributions with 
zero chemical potential and T^^ = 8 MeV, Tp^ = 5 MeV, and T^^ = 3.5 MeV, corresponding 
to average neutrino energies of {Ei,J = 25 MeV, {EpJ = 16 MeV, and {E^J = 11 MeV. 
More elaborate investigations of neutrino production in supernovae indicate that the high- 
energy tail of the neutrino spectra is better described by a Fermi-Dirac distribution with a 
finite chemical potential 

In what is considered the birth of neutrino astrophysics, neutrinos from supernova 
SN1987A have been detected by the Kamiokande and 1MB |^ water Cerenkov detec- 
tors (11 and 8 events, respectively). It is generally assumed that these events originated 
from the Uf. + p ~^ n + reaction in water. The detection of u^, and neutrinos via the 
u + e z/' + e' scattering or the ^^0(z/e, e~)^^F reaction was strongly suppressed by the small 
effective cross sections of these processes. The observability of supernova neutrinos will sig- 
nificantly improve when the Superkamiokande (SK) detector becomes operational [Q. This 
detector, with about 15 times the target material for supernova neutrinos of Kamiokande 
and a lower threshold of E^^ = 5 MeV, will be capable to detect also the recoil electrons 
from u + e u' + e' . In principle, induced neutrino-electron scattering events can be 
separated from everything else in SK using their angular distributions and energy spectra 
P). However, only about one third of the Ux + g scattering events will have energies distinctly 
larger than the recoil electrons from z^e + e and z/g + e scattering. Moreover, these higher 
energy electron recoils have to be separated by their direction from the much more numerous 
positrons from Uf. + p —>■ n + with the same energy. 

In this Letter we suggest another signal in water Cerenkov detectors which allows one to 
unambigously identify induced events. The basis of our proposal is the fact that SK can 
observe photons with energies larger than 5 MeV Schematically our detection scheme 
works as follows (Fig. 1). Supernova neutrinos, with average energies of ~ 25 MeV, 
will predominantly excite 1~ and 2~ giant resonances in ^^O via the ^^0{i>x, u'^Y^O* neutral 
current reaction ||T^. These resonances are above the particle thresholds and will mainly 
decay by proton and neutron emission. (Decay into the a channel, although energetically 
allowed, is strongly suppressed by isospin conservation |]10[)- Although the proton and 
neutron decays will be mainly to the ground states of ^^N and ^^O, respectively, some of 
these decays will go to excited states in these nuclei. If these excited states are below the 
particle thresholds in ^^N {E* < 10.2 MeV) or ^^O (E* < 7.3 MeV), they will decay by 7 
emission. As the first excited states in both nuclei {E* = 5.27 MeV in ^^N and E* = 5.18 
MeV in ^^O) are at energies larger than the SK detection threshold, all of the excited states 
in ^^N and ^^O below the respective particle thresholds will emit photons which can be 
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observed in SK. 

Of course, it is important to estimate the effective ^^0(z/^, z/^P7) and ^^O^u^, i^'^nj) cross 
sections for SK and to compare it to the effective "background" cross sections, stemming 
from the u^+p ^ n + e~^ and u + e ^ u' + e' events with energy release similar to the energy 
of the photons. Assuming equal luminosities for all neutrino species [i = i/g, J^e, ^x) leaving 
a supernova, the relative event rate (cxefr) for a specific neutrino-induced process in a water 
Cerenkov detector is 

n f 

(aeff) ~ 7T;r / dEf{E)a{E) (1) 



{E) 

where / is the neutrino energy spectrum and the factor 1/ {E) accounts for the ratio of fluxes 
for the different neutrino flavors. <j{E) is the total cross section for the neutrino-induced 
process and n is the number of targets for an individual neutrino process in a single water 
molecule (n = 10 for neutrino-electron scattering, n = 2 for z/g + p — n + e"*" and n = 1 for 
neutrino reactions on ^^O). 

To calculate the ^^0{h'x,J^'xPl) and ^^O^Uxj^'x^l) cross sections we assume a two-step 
process. In the first step we calculate the ^^©(i/^^, ^'^)^^0* cross section as a function of 
excitation energy in ^^O within the continuum random phase approximation (CRPA). In the 
second step we calculate for each final state with well-defined energy, angular momentum and 
parity the branching ratios into the various decay channels using the statistical model code 



SMOKER ||TT[, considering proton, neutron, a and 7 emission. As possible final states in the 
residual nucleus the SMOKER code considers the experimentally known levels supplemented 
at higher energies by an appropriate level density formula [Tl|]. If the decay leads to an 



excited level of the residual nucleus (e.g. to p-|-^^N*), we calculate the branching ratios for 
the decay of this level in an analogous way. Keeping track of the energies of the ejected 
particles and photons during the cascade, and weighting them with appropriate branching 
ratios and the corresponding differential ^^0(^'x, z/^)^^0* cross section, we determine the 
various particle (proton, neutron, a) and photon spectra for the (z/^, z/^) reaction on ^^O. 
We performed a similar calculation also for the [ux, v'^) reaction on ^^O. The contribution of 
various neutrino energies was weighted according to a (normalized) distribution f{E). Note 
that the same CRPA approach has been successfully applied to the muon capture on ^^O 



T4 . The model is described in details in Refs. ||13| , p!0| . As residual interaction we adopt the 
finite-range force based on the Bonn potential [Q. A similar two-step approach (combining 
CRPA and statistical model) has been tested successfully against the integrated (7^) / (777,) 
data on I'^O |T0|. 

The total and partial cross sections for Ux and Ux induced neutral current reactions on ^^O 
were evaluated first using the Fermi-Dirac neutrino spectrum with zero chemical potential 
and temperature T = 8 MeV (FDl) The results are listed in Table I. The total [ux, v'x) 
and [vxi ^'x) cross sections are roughly the same as the vector-axial vector interference term 
is rather unimportant. As expected, the partial cross sections for decay into proton and 
neutron channels dominate the total cross section. The proton channel is favored over the 
neutron channel by the lower threshold in ^^O. We find that a significant fraction of these 
decays go to excited states in ^^N and ^^O below particle thresholds and thus decay by 7 
emission. The relatively larger importance of this decay mode in ^^N {k, 24%) than in ^^O 
(~ 6%) reflects the larger number of final states in ^^N due to the higher particle threshold. 
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We find 3.2 ■ 10~^^ cm^ for tlie total 7 producing cross section for each flavor of Vx plus Vx 
in the neutral current reactions on ^^O, obtained by adding the {jrf) and (^7) partial cross 
sections. 

As discussed above, the Vx and Vx induced reactions will produce 7 events in the energy 
range ii^ ^ 5 — 10 MeV. The other events at these energies (a "background" for our purpose) 
will stem mainly from the z/g + p ^ + e"*" reaction. Adopting the Fermi-Dirac distribution 
with T = 5 MeV and zero chemical potential, we calculate a total cross section for this 
reaction of 47 ■ 10~^^ cm^. (This includes the factor n = 2 for the two protons in a water 
molecule. The result is somewhat smaller than that quoted in |^ where minor effects, such 
as the weak magnetism and recoil were not included.) However, the energy spectrum of 
positrons as seen by SK is peaked at around 15 MeV and only a small fraction of events are 
in the energy window E = 5 — 10 MeV. This becomes obvious in Fig. 2, where we compare 
the positron spectrum with the 7 spectrum calculated for the Ux and Ux induced reaction 
on ^^O. The latter has been multiplied by a factor 2 (to account for z/^ and u^- neutrinos) 
and by 16/25 to consider the ratio of Ue and z/^ fluxes (~ {E^J) / {EpJ, see Eq. (1))) at the 
detector. An energy resolution 14%/(E/10)i/2 ^^ere E is in MeV, i.e., 1 MeV for the 
energies of interest, has been assumed for the detector. As is obvious from Fig. 2, the 7 
spectrum constitutes a clear signal at E = 5 — 7 MeV on top of a smooth background from 
the i?e + p —>■ n + reaction. Our calculation predicts most of the photons to stem from the 
decay of the 3 lowest levels in ^^N and ^^O. Further, in Fig. 2 we assume that the detector 
will record all photons in a possible cascade of several 7 rays as a single event. Let us stress 
that each of such multi-photon events will contain at least one photon above the 5 MeV 
threshold. 

Other possible backgrounds are neutrino-electron scattering and charged current reac- 
tions on ^^O. For these reactions we find smooth electron or positron spectra, whose cross 
sections in the interval E = 5 — 10 MeV (normalized with the appropriate flux ratios and 
target numbers n) are much smaller than the 7 signal. Water also contains a tiny amount 
of ^'^O and even less ^^O. However, their natural abundances (~ 0.2% and 0.04%, respec- 
tively) are too small for neutrino reactions on ^^O to be of importance (see Ref. for the 
calculated charged current cross sections). 

We then repeated our calculation of the Ux and z/^ induced reactions on ^^O, using a 
Fermi-Dirac neutrino spectrum (FD2) with T = 6.26 MeV and /i = 3T 

~ TT ^ ' (2) 

1 + ex-p[{E - fi)/T\ 

which has the same average neutrino energy as the FDl distribution. We find that the FD2 
total and partial cross sections are smaller by about a factor of two when compared to the 
FDl results (see Table 1). Noting that the main contribution to the cross sections comes 
from neutrinos with E,y > {EJ)^ this scaling simply reflects the ratio of the two Fermi-Dirac 
distributions in that energy region. At the same time, the cross section for the dominant 
reaction z/g + p — > ^ + e"*" is changed only slightly, to 44 ■ 10~^^ cm^ when the FD spectrum 
with T = 5 MeV and /x = is replaced by a spectrum with T = 4 MeV and /x = 3T. 

We note that photodissociation data confirm a significant decay rate of the giant dipole 
resonance in ^^O by proton and neutron emission into excited states of ^^N and ^^O |T^. In 



agreement with our model these decays mainly lead to the first 3 excited levels in these nuclei 
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and are relatively larger in ^^N than in ^^O [0,^. While the total decay rate appears to 



be in reasonable agreement with our calculation, the data suggest a preference of the decay 
to the 3/2~ state at ~ 6.3 MeV over the decay to the positive-parity states at around 5.3 



MeV, caused by nuclear structure effects beyond our present model |[18|J19|| . This suggests 
that a fraction of the signal, predicted by our calculation at ~ 5.3 MeV just above the 
SK detection threshold, is to be shifted to 6.3 MeV, where it can be easier detected in SK. 
To clarify this point, a detailed investigation of the photodissociation process on ^^O within 
the current model is in progress. 

Superkamiokande is expected to detect about 4000 positrons from the Ue + P —>■ n + e'^ 
reaction |^ for a supernova going off at 10 kpc (~ 3 ■ 10^ lightyears or the distance to the 
galactic center). By scaling the respective effective cross sections, we estimate that such 
a supernova will produce about 360 (FDl) or 190 (FD2) 7 events in the energy window 
£" = 5 — 10 MeV, to be compared with a smooth background of about 270 positron events 
from the i^e+P n+e~^ reaction in the same energy window. This number of events produced 
by supernova neutrinos via the scheme proposed here is larger than the total number of 
events expected from i/a;-electron scattering (about 80 events pj). More importantly, the 7 
signal can be unambiguously identified from the observed spectrum in the SK detector, in 
contrast to the more difficult identification from i/j.-electron scattering. 

In conclusion, we propose a novel signal for the identification of supernova fi and r 
neutrinos in water Cerenkov detectors. Our suggestion is based on the fact that the levels in 
^^O, which are excited by inelastic neutral current scattering of supernova neutrinos, have 
decay branches via proton and neutron emission into excited states in ^^N and ^^O. These 
states, in turn, decay by emission of photons with > 5 MeV, which can be detected in 
the Superkamiokande water Cerenkov detector. We show that the expected number of 7 
events in the relevant energy window E = 5 — 10 MeV is noticeably larger than the positron 
or electron background expected from other neutrino reactions in water. It is amusing to 
note that the neutrinos from SN1987A at 50 kpc would have created about 10 photons 
in SK (which did not exist at that time, unfortunately), the same number of events as all 
of the recorded events in the Kamiokande or 1MB detectors then, which launched the era of 
neutrino astronomy. 

We are grateful to F.-K. Thielemann for providing us with his statistical model code 
SMOKER. We thank Y.-Z. Qian and F.-K. Thielemann for helpful discussions. This work 
was supported by the U.S. National Science Foundation (grants PHY94-12818 and PHY94- 
20470), by the U.S. Department of Energy (DE-FG03-88ER-40397) and by the Swiss Na- 
tionalfonds. 
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TABLES 



TABLE I. Total and partial cross sections for and induced reactions on °0, calculated 
for a Fermi-Dirac neutrino spectrum with temperature and chemical potential (T=8 MeV, /x = 0) 
(upper part) and (r=6.26 MeV, fx = 3T) (lower part). 



reaction 




reaction 




i60(i.„i.i)X 


5.90 


i60(P„i>i)X 


4.48 


160(z.„i/^p)15n 


3.75 


i60(P„ i>;p)i5N 


2.93 


I60(i/„i/»i50 


1.76 


i60(P„i>;n)i50 


1.29 


160(i/„i.;p7)15N 


1.41 


i60(i7„i7;p7)''N 


1.09 


i60(i/„i/>7)i50 


0.37 


i60(P„i>>7)i50 


0.28 




3.08 


i60(P„i>OX 


2.50 


160(z.„i/;p)15N 


2.02 


i60(i^„i>;p)i5N 


1.69 


I60(i/„z/»i50 


0.90 


i60(i^„P»i50 


0.70 


i60(i/„t'^P7)i5N 


0.72 


i60(i^„i^^P7)i5N 


0.59 




0.18 


160(77,. /y^j-oi'^'O 


0.14 
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FIGURES 



FIG. 1. Schematic illustration of the detection scheme for supernova and Uj- neutrinos in 
water Cerenkov detectors. 

FIG. 2. Signal expected from supernova neutrinos in a water Cerenkov detector. The solid 
line is the sum of the 7 spectrum, generated by and ij^ reactions on ^^O, and of the positron 
spectrum (dashed line) from the v'e+p ^ n + e'^ reaction. The upper part (a) has been calculated 
assuming Fermi-Dirac neutrino distributions with (T = 8 MeV, M = 0) and (T = 5 MeV ,^ = 0) 
for Vx and Ue neutrinos, respectively. In the lower part (b) Fermi-Dirac neutrino distributions with 
(T = 6.26 MeV, fi = 3T) and (T = 4 MeV, n = 3T) have been assumed for Ux and Ug neutrinos. 
The energy E refers to the photon or positron energy, respectively. The spectra are in arbitrary 
units. 
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